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In this paper we revisit an idea originally proposed by Mandelbrot about the possibility to observe 
"negative dimensions" in random multifractals. For that purpose, we define a new way to study 
scaling where the observation scale r and the total sample length L are respectively going to zero 
and to infinity. This "mixed" asymptotic regime is parametrized by an exponent x that corresponds 
to Mandelbrot "supersampling exponent". In order to study the scaling exponents in the mixed 
regime, we use a formalism introduced in the context of the physics of disordered systems relying 
upon traveling wave solutions of some non-linear iteration equation. Within our approach, we show 
that for random multiplicative cascade models, the parameter x can be interpreted as a negative 
dimension and, as anticipated by Mandelbrot, allows one to uncover the "hidden" negative part 
of the singularity spectrum, corresponding to "latent" singularities. We illustrate our purpose on 
synthetic cascade models. When applied to turbulence data, this formalism allows us to distinguish 
two popular phenomenological models of dissipation intermittency: We show that the mixed scaling 
exponents agree with a log-normal model and not with log-Poisson statistics. 



PACS numbers: Valid PACS appear here 

Multifractal processes are random functions (or mea- 
sures) that possess non trivial scaling properties. They 
are widely used models in many areas of applied and fun- 
damental fields. Well known examples are turbulence, 
internet traffic, rainfall distributions or finance. For the 
sake of simplicity we will consider only non-decreasing 
multifractal processes (which increments define a mul- 
tifractal measure) denoted hereafter M{t). In the se- 
quel M{I) will stand for the measure of the interval /, 
M (/) = Jj dM and M(i, £) =M{[t,t + i]). Multifractals 
are characterized by the scaling of the partition functions: 
If one covers the overall sample interval of length L with 
L/£ disjoint intervals of size i, {Ii{i)}i=i...L/e: one usu- 
ally defines the order q partition function which scaling 
behavior defines the exponent TQ{q): 



L/e 



Z{q,i) = yM[lS)]' - r°(«) 



=1 



£-^0 



(1) 



where the limit ^ — > simply means that i/T —> 0, T 
being the correlation scale usually referred to as the in- 
tegral scale. When ro(g) is a (concave) nonlinear func- 
tion of q, the measure M{t) is said to be multifractal 
or intermittent. Within the multifractal formalism in- 
troduced by Parisi and Frisch (see e.g. [l|), the non- 
linearity of To{q) is interpreted in terms of fluctuations of 
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pointwise singularity exponents of the measure. Indeed, 
according to this formalism, tq (q) is obtained as the Leg- 
endre transform of the singularity spectrum fo{a) that 
gives the (Haussdorf ) dimension of the sets of points t 
of singularity a {M{t,i) ~ £"). Therefore q can be in- 
terpreted as a value of the derivative of fo{a) and con- 
versely a is a value of the derivative of To{q). Multi- 
fractality is also closely related to the notion of stochas- 
tic self-similarity: The measure M(t) is self-similar in a 
stochastic sense if, for all s < 1, M{st)=iaw sWsM{t) 
where Wg is a positive random weight independent of 
M. It can be easily shown that this stochastic equality 
implies that the expected value of Z{q,£) (i.e., the order 
q moment of the measure) behaves as a power law, i.e., 
E[Z{qJ)] = Cqr(«) where T{q) is nothing but the cu- 
niulant generating function of In W^ ■ If one uses the ter- 
minology introduced in physics of disorderded systems, 
the exponent T{q) ~ InE [Z((7,£)] is defined from an "an- 
nealed" averaging while To{q) ~ E [lnZ((7,£)] is the ana- 
log of a free energy computed as a "quenched" average. 
As it will be discussed below, these two functions can be 
different for large values of q. 

The paradigm of self-similar measures are random mul- 
tiplicative cascades (for the sake of simplicity we will ex- 
clusively focus, in this paper, on discrete cascades but all 
our results can be easily extended to recent continuous 
cascade constructions ^j S H) originally introduced by 
the russian school for modelling the energy cascade in 
fully developed turbulence and to which a lot of math- 
ematical studies have been devoted [a, H S Q ■ Let us 



summary the main properties of these constructions. The 
integral scale from which the cascading process "starts" 
is denoted as T. A dyadic discrete cascade is build as 
follows: A measure is uniforniely spread over the start- 
ing interval [0, T] and one splits this interval in two equal 
parts: On each part, the density is multiplied by (posi- 
tive) i.i.d. random factors W such that E \W] — 1. Each 
of the two sub-intervals is again cut in two equal parts 
and the process is repeated infinitely. It is convenient to 
introduce the random variable w = In(W^) which prob- 
ability density function will be denoted as g{x). Thus 
g{x) is Gaussian or Poisson for respectiveley log-Normal 
and log-Poisson cascades. The stochastic self-similarity 
property of the limit measure M associated with pre- 
vious iterative construction can be directly proven and 
it is easy to show that r(q) is related to the cumulant 
generating function of uj: 



to define a new exponent T^{q) as follows: 



T(g)=g-ln2(E[e''"])-l 



(2) 



Moreover, as first established by Molchan [T, 9], the mul- 
tifractal formalism holds for random cascades and one 
has the following relationship between ro(g) and T(g), for 
g > 0, (note that a similar relationship holds for negative 
values of g [3]): 

To{q) = T{q) if q < go and Tn{q) = a^q if g > go (3) 

where go is the value of g corresponding to minimum 
value of /o(a) (in general /o(ao) = 0). This discrepency 
between the annealed and quenched spectra has been ex- 
tensively studied on a numerical ground in ref. [lO| and 
was referred to as the "linearization effect". It has no- 
tably been observed to be independent of the nature of 
the cascade and of the overall length L of the sample. As 
it will be emphasized below this effect has not been prop- 
erly taken into account in the literature of turbulence (see 
e.g. P). 

One of the goals of this paper is to recover the lin- 
earization effect and to establish how it can be somehow 
controlled. For that purpose, we will consider the scal- 
ing of partition function ([1]) in some "mixed" asymptotic 
regime where, as the resolution becomes smaller, the to- 
tal length of the sample is inscreased. Let us introduce 
some useful notations: We call I the scale of observations 
(it can be the sampling scale or a multiple of it), T the 
integral scale and L the total sample length. Ni will refer 
to the number of samples per integral scale and Nt the 
number of integral scales. We have obviously Ni — T£^^ 
and Nt = LT^^ and N ^ NiNt = L/£ is the total num- 
ber of samples. If T is fixed, the limit N — > -t-oo can 
be conveniently controlled using an additional exponent 
X > as Nt ^ Nf. Let us mention that such an ex- 
ponent has been already introduced by B.B. Mandelbrot 
as an "embedding dimension" [l^, [l3l in order to discuss 
the concept of negative dimension and latent singulari- 
ties (see below). Within this framework, the definition 
P]) of the partition function depends on x and allows us 



Lr 



Z(g,£,x)=5]M[/,(^)]'^ r-(«)' 



(4) 



The two "extreme" cases are : (i) the x = case which 
corresponds to a fixed number of integral scales Nt while 
I —f and (ii) the x = +oo case which corresponds 
to a fixed observation scale £ while Nt -^ +oo. It 
results that for x = one recovers former definition 
(UJ) of ro(g) while roo(g) = T(g) as defined in ([2]). In 
that respect x allows us to interpolate between quenched 
and annealed situations. In order to compute T-^{q) one 
needs to study the behavior of the probability law of 
Ziq,£,x) [l3]- For that purpose, along the same line as 
in references [la, [lali 1^* us study its Laplace transform: 



G{s,* 



E 



'Z(q,l,x) 



Let r,p two integers and 



let us define the iteration ?7i — > m + 1, i ^ 2~p£ and 
L — > 2^L. In other words, at each iteration step, the 
resolution is divided by 2^ while the number of indepen- 
dent integral scales is multiplied by 2''. One thus have 
Nt = 2™ and Ne = 2^"^ that corresponds to x = r/p. 
Within this parametrization, G{s,i) will be denoted as 
G{s,m,p,r) and G{s,m, 1,0) will be denoted as H{s,m). 
If Af is a random cascade as defined previously, its self- 
similarity allows one to prove that G{s,m,p,r) can be 
written as |17| : 



G{s,m,p,r) = {H{s,pm)) 
where H(s,m) satisfies the following recursion 



(5) 



H{s,m+l)^[H{s,m)*g{s + lii2)Y (6) 

where g{x) is the pdf of lo the logarithm of cascade 
weights and * stands for the convolution product. 

It is easy to see that Eq. ([6|) as two uniform "sta- 
tionnary' solutions G{s,n) — and G{s,n) = 1, the 
first one being stable while the latter is (linearly) un- 
stable. The initial condition connects the stable state 
H{—oo,0) = to the unstable one H{oo,0) = 1 and, 
as shown in [l^ [l^, [l^ , this kind of non- linear equation 
admits traveling wave solutions H{s,n) = Hq{s — vn) 
where v is the front velocity. This velocity can be com- 
puted by studying the linearized version of Eqs. ([5]) 
and ([6j) around the unstable state. After a little al- 
gebra [12], one can show that the solutions of ([5|) are 
traveling fronts that can be written, when s — > -f oo, as 
G{s,m,p,r) = (1 - C'e-T(^-''(^)P"))2™, provided ^(7) 
satisfies the dispersion relation: w(7) = "t ){r-pr{'r)) 
where t(7) is defined by Eq. ([2]). One can reproduce 
the same kind of analysis as in refs. [l3, 120, |2l| and 
show that a standard Aronson- Weinberger stability cri- 
terium can be used to compute the selected velocity and 
7 values: Let g^ be the unique positive 7 value such that 
v{q^) = min^>o "^(7). The selected velocity 11(7) actually 
corresponds to the 7 value equal to the exponential de- 
creasing rate of the initial condition provided it is greater 
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FIG. 1: Tx(q) estimated for log-normal random cascade for 

X = (o), 0.2 (x), 0.5 (A), 0.7 (D), 1 (0), 1.5 (v). For 

each value of x, we have chosen T = 4096, A^ = 0.2 and 

Nt depends on x such that NrNg is fixed. One clearly sees 

that, as the value of x increases, the value of q^ increases, the 

slope a^ of T^{q) decreases and T-^{q) becomes closer to the 

annealed spetrum r(g) (solid line). The dashed lines represent 

the analytical expressions of r^ derived from Eqs. ^ and Q . 

In the inset the same estimations have been performed for a 

log-Poisson cascade. As expected, in that case Tx{q) weakly 

depends on x and rapidly converges towards r(g) represented 

by the solid line (see text) 



than q^. But the initial condition iJ(s,0) is precisely 
given by the Laplace transform of the unconditional law 
of M[0, 1]9 and if E [M«] = Mq < +00 (such a condition 
is well known to be satisfied provided T{q) > 0) then 
H{s, 0) ~s^oo 1 — MqC''^". Therefore, the selected veloc- 
ity is simply v — - )i^-P'^w)\ if g < g^ and v = v{q^) 
otherwise. Thanks to the fact that m is related to the 
resolution scale by ln(£) — —mp\ii(2), and using x = "^ IVi 
the velocity, measured as respect to In(^), finally becomes 
Iktx ,-f „ / „ ^^A .,(r,\ - max„>n ^^^2^ other- 



^(9) = „ if 9 < "Zx ^'^'^ ^(9) 



(j>0 



wise. If one denotes /(a) the legendre transform of r((7), 
the value q-^ for wich the maximum value of '^^'^'^^ is 
reached, satisfies q^ = f'{a^) with 

/("x) = -Xandax =T'(qx) =w(gx) (7) 

Since positive values of /(a) correspond to a fractal di- 
mension, the X value can be seen as a kind of "nega- 
tive dimension" . In order to solve the initial problem 
one can reproduce the analysis of [l^l to show that the 
front velocity is directly related to the mean value of 
\nZ{qJ,x) as q-^K[\a.Z{q,i,x)] ~ v{q)hi{t). Since 
in the "moving frame" , the probability distribution of 
\nZ{q,i,x) converges, it results that the fluctuations 
of ln2'(g, £,x)/ln(€) vanish. This notably implies that 
q-i " (g, ,x) _^ v{q), where the convergence is in prob- 
ability. According to definition (|4]), this is equivalent to 
say that T^{q) is non random and equals to v{q,x) + Xi 
i.e.: 



-I- 




'^xil) = '^(9) if 9 < 9x 
T^{q) — qa^ otherwise 



(8) 
(9) 
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FIG. 2: Same plot as in Fig. [T] but for the dissipation field 

e{x) estimated from a high Reynolds number jet turbulence 
experiment. The symbols have the same meaning as in Fig. [T] 
As for the log-Normal synthetic cascade one clearly observes, 
as X increases, a continuous evolution of Tx{q) towards the an- 
alytical parabola corrresponding to a log-Normal model with 
A'^ = 0.2 (solid line). Such an effet is not observed for a 
log-Poisson cascade. 



where a^ is defined in ([7]). Let us note that a rigourous 
proof of Eqs ([8]) and ([9]) will be provided in flil] . In the 
case X = o^^ recovers standard linearization effect (Eq. 
dSl)) which has been generalized to any value of x in the 
mixed asymptotic regime. As x increases so does q^ and 
T^(g) continuously converges towards T{q). Singularities 
a < ao have been qualified by Mandelbrot as "latent" 
because they are only observable for large enough values 
of the " supersampling" exponent x II4I • 

In order to illustrate our results, we have performed 
several numerical simulations on both continuous and 
discrete cascades which statistics are log-normal (LN) 
and log-Poisson (LP). In the log- normal case one has 
T{q) = q{l + A^/2) - A2g2/2 - 1 and therefore by solv- 
ing Eq dZl) one gets a^ = 1 + X'^/2 - \^2{+x) and 

g^ = X~^y/2{l + x)- If Fig. [U we have plotted, for 
each value of x = 0, 0.2, 0.5, 0.7, 1, 1.5, C^{q) = 1 + T^{q) 
as a function of q estimated from the scaling of par- 
tition function Z{q^ ^, x) over a wide enough range of 
scales. The log-normal cascade samples are generated 
from discrete cascades with T = 4096, A^ — 0.2 and, 
for each value of x^ ^t is chosen such that, at the 
smallest scale, the number of sampling points is almost 
constant. As expected, one sees that, when the value 
of X increases, q^^ increases, a^ (the slope of T^(g) for 
q > q^) decreases and one has Ty^{q) = T{q) over a wider 
range of q. The same kind of simulations have been 
performed on log-Poisson synthetic cascades for which 
T{q) ^ q{l + \^{e^ ^ l)/5^) + \^{1 - eiyS^) - 1. In that 
case, it is easy to show that since T{q) has an asymp- 
tote when q -^ -fcx) t' (q) > 1 + \'^{e^ - 1)/S'^) and 
/(a) > ~X^/S^. Therefore, for all x > ^ 1^^^ one has 
'^x^'i) ~ '^{l)- ff one chooses the log-Poisson parameter 
values usually considered to model energy dissipation in 
turbulence [i|,[24|i i-^-i ^^ = 0-2 and 5 = ln(2/3), one sees 
that T^{q) rapidly converges towards r(q) and no longer 
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FIG. 3: Plot of a^, the slope of r^iq), for q > Qx estimated 
for the turbulent dissipation field as a function of ^1 + x- 
The straight solid line represents the log-Normal prediction 
while the dashed curve is the log-Poisson prediction. Error 
bars have been estimated using Monte-Carlo trials of the Log- 
normal cascade. This representation allows us to distinguish 
clearly the two models of energy dissipation intermittency. 



varies for x > 0.2. This is illustrated in the inset of Fig. [T] 
where all the estimated T^{q) for x = Oi 0-2, 0.5, 0.7, 1, 1.5 
are close or equal to the annealed spectrum T{q) repre- 
sented by the solid line. 

According to our analysis and from numerical experi- 
ments reported in Fig. [U one sees that despite the fact 
that T-)^ (q) associated with log- normal and log-Poisson are 
very similar for small values of q and small values of x, 
the situation is very different for large values of x ^nd 
q. Since both models are traditionnally used to describe 
the spatial flucuations of energy dissipation in fully de- 
veloped turbulence [l|, \2^ , we naturally reproduced pre- 
vious statistical analysis using experimental data of tur- 
bulence. The data have been recorded by the group of B. 



Castaing in Grenoble in a low temperature gazeous He- 
lium jet experiment which Taylor scale based Reynolds 
number is R\ = 929 [23] . Assuming the validity of Tay- 
lor hypothesis and isotropy of the flow, a proxy of the 
dissipation field e{x) ~ (dv/dt)^ is build from the tem- 
poral longitudinal velocity signal. The overall sample is 
such that Nt — 2500. The estimated values of T^{q) are 
reported in Fig. [21 The similarity with the results of 
log-normal synthetic cascades is striking: As in Fig. [T] 
when X increases, the estimated T^{q) converges to the 
paraboHc T{q). In Fig. [3l we have plotted the values 
of the asymptotic slope of T^{q), a^^ ^s a function of 
^1 -f X'- one clearly sees that the data perfectly match 
the linear log- Normal expression (solid line) which is very 
different from the log-Poisson prediction (dashed line). 

To summarize, we have shown in this paper, using trav- 
eling wave solutions of cascade non-linear iteration equa- 
tion, that quenched and annealed averaged partitions 
function have different behavior for values of q larger 
than a critical value q^, analog of the glass transition 
temperature in spin glass systems. This difference can be 
controlled using some " supersampling" exponent x which 
defines an asymptotic limit that mixes small scales and 
large number of samples regimes. By analyzing the value 
of the multifractal scaling exponents for various values 
of X 0116 can distinguish between different cascade mod- 
els. In the context of the modelling of energy dissipation 
intermittency in fully developed turbulence, we have pro- 
vided evidences supporting log-Normal statistics against 
log-Poisson statistics. 
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